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ABSTRACT 

We use the redshift distribution of type-la supernovae (SNe) discovered by the 
Supernova Cosmology Project to constrain the star formation history (SFH) of the 
Universe and SN la progenitor models. Given some of the recent determinations of the 
SFH, the observed SN la redshift distribution indicates a long lh~^ Gyr) mean 
delay time between the formation of a stellar population and the explosion of some of 
its members as SNe la. For example, if the Madau et al. (1998) SFH is assumed, the 
delay time r is constrained to be r > 1.7(r > Q.7)h~^ Gyr at the 95%(99%) confidence 
level (CL). SFHs that rise at high redshift, similar to those advocated by Lanzetta et 
al. (2002), are inconsistent with the data at the 95% CL unless r > 2.5h~^ Gyr. Long 
time delays disfavor progenitor models such as edge-lit detonation of a white dwarf 
accreting from a giant donor, and the carbon core ignition of a white dwarf passing 
the Chandrasekhar mass due to accretion from a subgiant. The SN la delay may be 
shorter, thereby relaxing some of these constraints, if the field star formation rate falls, 
between z = 1 and the present, less sharply than implied, e.g., by the original Madau 
plot. We show that the discovery of larger samples of high- 2; SNe la by forthcoming 
observational projects should yield strong constraints on the progenitor models and 
the SFH. In a companion paper, we demonstrate that if SNe la produce most of the 
iron in galaxy clusters, and the stars in clusters formed at z ~ 2, the SN la delay time 
must be lower than 2 Gyr. If so, then the Lanzetta et al. (2002) SFH will be ruled out 
by the data presented here. 
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1 INTRODUCTION 

The evolution of the SN rate over cosmic time is a key in- 
gredient for understanding the chemical enrichment history 
of the Universe and the formation and properties of galaxies 
and clusters. Since the progenitors of core-collapse SNe are 
massive, short-lived stars, their cosmic rate evolution will 
closely follow the universal star formation history (SFH). On 
the other hand, the actual route that leads a white dwarf 
(WD) to explode as a SN la is still an open question (see, 
e.g., Yungelson & Livio 2000, and references therein). 

The mechanism that leads to the explosions of SNe la 
likely involves a significant delay between the formation of 
the progenitor system and the explosion of the SN. The rate 
of SNe la at a given epoch is therefore expected to depend 
on the star formation rate up to several Gyr preceding that 



epoch. Thus, the cosmic evolution of SN rates can be used 
to constrain both the global SFH, and the characteristic de- 
lay time between star formation and SN la explosion (e.g., 
Madau, DeUa VaUe & Panagia 1998; hereafter MDP). Since 
competing SN la models predict different delay times, obser- 
vational constraints on the delay can discriminate between 
some of the proposed scenarios for SN la progenitor systems. 

During the last few years, SNe la have been used suc- 
cessfully as cosmological distance indicators (e.g., Riess et 
al. 1998,2001; Perlmutter et al. 1999; Tonry et al. 2003). Per- 
haps the main uncertainty that still plagues SN-Ia-based dis- 
tance measurements is the possibility of evolution in SN la 
properties. The use of SNe la as distance estimators relies 
on the assumption that distant SNe la at redshifts as high 
as 2: = 1.7 are similar to local events. In order to study this 
question from a theoretical perspective, we require knowl- 
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edge on the nature of the progenitor systems of SNe la. For 
any given scenario, it is possible to estimate the evolution 
of high-^; systems relative to local ones, e.g., as a result of 
metallicity or stellar-age effects. Determination of SN la de- 
lay times, with the resulting constraints on progenitor mod- 
els, are therefore desirable. 

Motivated by these questions, a number of authors have 
recently studied the evolution of SN rates. Most studies have 
combined models of the SFH with a recipe for the delay 
function of SNe la, and have calculated the expected evo- 
lution of SN rates, either per unit comoving volume, or per 
unit stellar luminosity (e.g. J0rgensen et al. 1997; Sadat et 
al. 1998; MDP; Ruiz-Lapuontc & Canal 1998; Yungclson & 
Livio 1998; Kolatt & Bartelmann 1998). As measurements 
of SN la rates at high z have became available (Pain et al. 
1997; 2002) attempts have been made to compare the ob- 
served rates with model predictions (e.g., MDP; Sadat et 
al. 1998; Kobayashi ot al. 1998; Kobayashi, Tsujimoto, & 
Nomoto 2000; Pain ct al. 2002; Calura & Matteucci 2003). 
However, no strong conclusions have been reached regarding 
the SFH or SN la progenitors. 

All available measurements of SN rates beyond the local 
Universe (Pain ot al. 1997; 2002; Hardin ct al. 2000; Gal- 
Yam, Maoz, & Sharon 2002; Tonry et al. 2003) are based 
on SNe la in the redshift range z = 0—1. Due to the 
limited number of observed SNe la, these studies have dis- 
tributed their samples into wide redshift bins. In Gal- Yam 
ct al. (2002) we divided our sample into two redshift bins, 
with < z >= 0.25 and < z >= 0.9. This binning caused no 
loss of information in our work, as the number of bins was 
comparable to the number of SNe (2 — 3). On the other hand. 
Pain et al. (1997; 2002) and Tonry et al. (2003) used only 
one bin, and calculated the SN rate at an average redshift, 
z ~ 0.5. In these cases, most of the redshift information is 
lost in the binning process. 

This loss may bo averted if, instead of deriving absolute 
rates from the observations, and comparing them to pre- 
dicted rates, we begin with the predictions and fold them 
through the observational filters. We can then compare the 
prediction for a particular experiment with the unbinned 
observations. Specifically, wo can test whether or not the 
unbinned distribution of SN redshifts in a particular sur- 
vey is consistent with some model. Indeed, several authors 
have calculated such distributions (e.g., Dahlen & Prans- 
son 1999, hereafter DF; Sullivan et al. 2000a). DF noted 
that these distributions can be used to constrain progenitor 
models. However, they concentrated on future, deep surveys 
with the James Webb Space Telescope, and did not attempt 
to constrain model parameters using existing SN data. 

Pain et al. (1997, 2002) and Tonry et al. (2003) have cal- 
culated the expected redshift distributions of SNe in their 
respective surveys, and have compared them with the ob- 
served distributions Eis part of their derivation of SN rates. 
However, the SN redshift distributions calculated by Pain 
et al. and Tonry et al. are based on particular assumptions 
about the evolution of the SN la rate between z — and 
z ^ 1, which is assumed to be either constant (Tonry et 
al. 2003), or to vary as some power of the redshift (Pain 
ct al. 2002). The comparison between the expected and ob- 
served SN la redshift distributions is used to find the best 
fitting rate evolution, and to derive the average SN la rate 
at the mean redshift {z ~ 0.5 in both cases). Again, the in- 



formation contained in the redshift distribution of SNe la at 
z = — 1 is lost in the averaging process. Both groups note 
that their data seem to indicate a slowly varying SN la rate 
in this redshift range. Tonry et al. conclude that, given that 
the slope of the SN la rate is shallower than that of some 
determinations of the SFH in this redshift range, the typical 
delay time between star formation and SN la explosion must 
be ~ 1 Gyr. 

In view of the fact that some models (e.g., MDP) predict 
strong evolution in the SN la rate at 2: = — 1, in contrast 
with the trends seen by Pain et al. (2002) and Tonry et al. 
(2003), we have undertaken a more comprehensive investi- 
gation of the SN la redshift distribution in this range. As 
we will show below, existing data can already place interest- 
ing limits on the SFH and on the SN la characteristic delay 
time. 

Throughout this paper, we assume a flat cosmology with 
fim = 0.3 and Qa = 0.7. We designate with h the Hubble 
parameter in units of 100 km s"^ Mpc"'^. 



2 CALCULATION OF THE SUPERNOVA 
REDSHIFT DISTRIBUTION 

We begin by calculating N{z, Rum), the redshift distribution 
of SNe visible in a region of sky at a given moment, to a 
limiting i?-band magnitude Rum- Throughout this paper, 
wo assume that all apparent magnitudes are measured in 
the i?-band. Our calculations can be applied to any other 
band, as long as all apparent magnitudes are measured in 
the same band. Then, 

N{z,Riim)^n{z)xdViz)xT{z,Rii,^) , (1) 

where n{z) is the SN rate density, defined as the number 
of SNe per unit time per unit volume as a function of the 
redshift, and dV{z) is the volume between z and z + dz 
in that region of sky. T{z, Rum) is the effective visibility 
time of a SN at redshift z, given the detection efficiency 
as a function of magnitude (this is often referred to as the 
"control time" ) . T depends on z because SNe at higher z are 
fainter, and therefore can be detected for a shorter period 
above the limiting magnitude Rum- The brightness of the 
SN in the chosen bandpass also depends on z since different 
parts of the SN spectrum are redshifted into the observed 
band for SNe at different redshifts. 

The SN la rate density, n{z), is a convolution of the 
star-formation history (SFH) with a "delay" or "transfer 
function" , which is the SN la rate vs time following a brief 
burst of star formation. The delay function accounts for the 
time span between star formation, through stellar evolution, 
the formation of WDs, and the final stage where a WD ac- 
cretes material from (or spirals in and merges with) a binary 
companion (MDP; DF). We follow here the delay function 
parameterization given by MDP. The overall time delay in- 
cludes the mass-dependent lifetime of the progenitor as a 
main-sequence stax, AIms- Once the progenitor has gone 
off the main-sequence and has become a WD, it has a prob- 
ability cx exp(— ■^) to explode as a SN la, where At is the 
time since the star left the main sequence. DF have pro- 
posed a simpler model, with a discrete delay time, t, which 
we further consider in section § 4.2. Since we will be in- 
terested solely in the redshift distribution of SNe, we will 
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ignore the normalization of n{z). Integrating the explosion 
probabilities over the stellar initial mass function, dN/dm, 
and the past SFH, ^'(t), the SN rate density as a function 
of cosmic time is 



t/ m, 



n{t) oc / *(t')dt'x 
Jo 

. t-e-AtMS dN 
exp( ) -; — (m) dm , (2) 

dm 

where, mmin and mmax are the minimum and mgiximum 
initial masses that will lead to the formation of a WD that 
explodes as a SN la. Following MDP, we adopt 

t — t' 

mmin = max[3M0, (— — ;q — )~°'^Mq], rrimax = 8Mq, 
iU vjyr 

AtMS I m 2.5 ir,\ 

^"'^ lOG^ = (mT) • (3) 

A Salpeter (1955) initial mass function (IMF), dNjdm ~ 
m-2■3^ is assumed. We convert n{t) to n{z) using the trans- 
formation between cosmic time and redshift. 



At = 



Ho 



dz 



{1 + Z)[(l+Zf^m + ^KV' 



(4) 



The detailed form of the SFH, ^, is a much debated 
issue in recent years. Initial studies of the UV to IR lumi- 
nosity density of star- forming galaxies (e.g., Lilly et al. 1996; 
Madau et al. 1996; Connolly et al. 1997, Madau, Pozzetti, & 
Dickinson, 1998) suggested that the SFH increases sharply 
between rcdshifts of zero and 2: ~ 1 — 2, and then decreases 
at higher z. However, this picture has been challenged. For 
instance, Steidel et al. (1999) argue for an almost-constant 
SFH sX z ^ 1 — A, while Lanzetta et. al. (2002) favor a 
monotonic increase in SFH out to z = 10. At lower rcdshifts, 
Cowio, Songaila, & Barger (1999) and Sullivan et al. (2000a) 
have found the SFH slope at « ~ — 1 to be quite shallow, 
^(z) ~ (1-1-^;), but recent measurements by Hippelein et al. 
(2003) seem to confirm the sharp rise, '9 ~ {l + z)"^, initially 
reported by Lilly ct al. (1996). 

To parameterise this range of possibilities for ^{z), we 
represent the SFH by a broken power law, with 
* {z) ~ (1 + 2)" at high z, and *(«) ~ (1 -I- z)'' at 
low z. These functions are joined smoothly at 2; = 1.2, using 
the prescription of Beuermann et al. (1999; see also Bersier 
et al. 2003), 



axs _|_ ( 



X ^(1.2) 



1 + z 



ll/s 



(5) 



with s = 5. The high-z index a assumes values between a = 

—2, which corresponds to the results of Madau ct al. (1998), 
and a = 2, which describes the work by Lanzetta et al. 
(2002). Similarly, the low- 2: index f) varies between /3 = 4, as 
advocated by Hippelein et al. (2003), and /3 = 1, tracing the 
results of Cowie et al. (1999)*. Fig. la shows two examples 
of the SFH parameterization. With these ingredients, each 
model for the SN rate density n(z) is determined by three 



* In the case of /3 = 1, the Beuermann prescription causes the 
high-2 portion of the SFH function to remain somewhat affected 
by the low-z Index /3 even at redshifts much larger than z = 1.2. 
In that case, we therefore use a simple broken power law. 



parameters - the typical delay time r and the two indices 
describing the SFH, a and !3. Fig. lb shows two examples 
of SN la delay functions, and Fig. Ic shows two examples of 
n{z) for particular combinations of SFH and r. 

For our chosen cosmology, the volume element is given 

by, 



dV{z) oc Da x c — oc 



dt [J^xdzf 



(l-h2)2 (I + Z) 



(6) 



where Da is the angular diameter distance, and x is given 
by 



X = 



[Qm{l + Z')S+QAV/^ 



(7) 



In the calculation of the effective visibility time 
T{z,Riim) we follow the procedures described in Gal- Yam 
et al. (2002). Briefly, assuming initially that SNe la are per- 
fect standard candles with identical light curves, the peak 
magnitude of a SN at a given redshift is Rpeak{z), and the 
period it spends above some limiting magnitude Rum can 
be easily determined. The efficiency of a survey, r], is gen- 
erally a function of the apparent magnitude R (with Rum 
defined by rj^Rum) = 0), and is assumed to be independent 
of the SN redshift, i.e., the detection probability is the same 
for SNe with the same magnitudes but different redshifts. 
Weighting the time the SN spends above the detection limit 
by the relevant efficiency values, gives 



T{z) 



-f 

Jo 



v[Rit, z)]dt 



(8) 



where R{t, z) is the light curve of a SN la at redshift z, 
and 77 = when R{t) > Ru,n- Note that cosmological time 
dilation will slow down the evolution of the light curves (e.g., 
Lcibundgut et al. 1996, Goldhaber ct al. 1997, Riess et al. 
1997), and hence lengthen the visibility times of SNe la by 
(1 + 2), but this is canceled out in the calculation of A'^(2) 
by the (1 + 2) reduction of the observed SN rate at redshift 
z. 

The SN survey data whose redshift distribution we will 
model were obtained in the R band, and range in redshift 
from 2 = 0.361 to 2 = 0.778. We require information on 
R{t) and Rpeak{z) in order to calculate T{z). To determine 
Rpeakiz), we fit a 3rd-order polynomial to the observed R- 
band peak magnitudes of high-2 SNe reported by Perlmutter 
et al. (1998, 1999). Note that no K-corrections are needed, 
as the input we require are the peak R magnitudes at all red- 
shifts, and these do not have to be translated into rest frame 
B magnitudes, as done by the SCP for Hubble-diagram pur- 
poses. In terms of the light curves, R{t), the observer-frame 
R band is a fair match for rest-frame B band in most of the 
redshift range we model. However, at the highest and low- 
est redshifts, rest-frame U and V wavelengths, respectively, 
dominate the observer-frame R band. Since the light curves 
of SNe la have slightly different forms in each band, we have 
constructed a set of 2-dependent observer-frame light curves, 
by linearly interpolating between the appropriate rest-frame 
light curves. The template light curves were adapted from 
U, B and V light curves, kindly supplied by B. Leibundgut, 
and supplemented with data points taken from Riess ct al. 
(1999), which also give slightly better coverage of the pre- 
maximum period. Early Z7-band observations of SNe la are 
scarce, so we had to extrapolate the rest-frame U -band light 
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redshift (z) redshift (z) 

Figure 1. Illustration of the modeling and comparison to data. Panel a shows two examples of the SFH, ^(2:) - a "Madau" SFH, with 
a peak at 2 = 1.2 (solid curve), and a shallower model (dashed curve) reflecting the proposed modifications by Cowic et al. (1999) and 
Steidcl ct al. (1999). Panel b shows two examples of the expected SN la rate density following a brief burst of star formation. These 
delay functions are calculated using the MDP prescriirtion, with characteristic exponential delay times of r = Ih^^ Gyr (solid) and 
T = Zh^^ Gyr (dashed). For display purposes, an arbitrary redshift of 2 = 3 has been chosen for the burst of star formation. SFH 
models are convolved with a delay function, and the resulting SN rate densities n{z) for a "Madau" SFH with r = Ih^^ Gyr (solid) 
and a "Cowie-Steidel" SFH with r = 3/i~^ Gyr (dashed) are shown in panel c. Panel d shows the predicted SN distributions, N(z), for 
the models of panel c, in a survey with the same observational parameters of the SCP search described in § 3. KS tests show that the 
cumulative version of N{z) from a model combining a "Madau" SFH with a typical delay time of r = lh~^ Gyr (panel e) is ruled out 
by the data, while a model with "Cowie-Steidel" SFH and r = Zh~^ Gyr is consistent with the data (panel f). Vertical axis units are 
arbitrary in panels a-d. 



curve. However, this extrapolated curve is only relevant for 
high-2 SNe, which are invariably discovered close to the hm- 

iting magnitude. Therefore, the influence of the extrapolated 
part of the hght curve, well below the peak magnitude, is 
negligible. The Rpeakiz) curve and observer frame _R-band 
light curves represent the properties of an average SN la. 

In reality, the peak magnitudes and light curve shapes 
are not uniform. "Peculiar" SNe la arc quite common in low- 
z SN samples, but they are apparently absent from high-z 
samples like the ones we consider (Li et al. 2001). One may 
worry that many SNe la of the "underluminous" , 1991bg- 
like, variety are missed at high z due to their lower luminos- 
ity, and we return to consider this possibility in section § 4.2. 
Even spectroscopically "normal" SNe la arc not perfect stan- 
dard candles, and exhibit an intrinsic scatter of 0.2 — 0.3 mag 
in peak luminosity (see Branch 1998, for a review). The peak 
magnitude of a SN la is correlated with the light-curve shape 
of the object, with brighter events having a slower rise to 
maximum followed by a slower decline (e.g., Phillips 1993; 
Phillips et al. 1999; Riess, Press, & Kirshner 1996; Riess et 



al. 1998; Perlmutter et al. 1995, 1999). A correction factor 
can be applied to the peak magnitude to calibrate SNe la as 

standard candles. The intrinsic scatter in peak magnitudes 
and the differences in light curve shapes may cause a scat- 
ter in the visibility times of SNe, since brighter and broader 
SNe will have longer visibility times. 

Perlmutter et al. (1999) show the distribution of correc- 
tion factors (their Fig. 4). In order to account for this distri- 
bution, we have calculated the effective visibility times for 
SNe with various peak magnitudes and light-curve shapes. 
We then average these times using the distribution of peak 
magnitudes given by Perlmutter et al. (1999) as a weighting 
function. We thus obtain a function T{z, Rum) that fits a 
population of SNe la with properties similar to those ob- 
served by the SCP. 

Having determined T{z,Riim) for a particular survey, 
we can calculate the redshift distribution of SNe la, N(z), 
as a function of the model parameters. Figure Id illustrates 
N{z) for the two examples of n{z) in Figure Ic. We then 
convert these redshift distributions to cumulative form, and 
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Table 1. The SCP SN sample 



SN 


Redshift 


1997cl 


0.636 


1997em 


0.460 


1997ep 


0.462 


1997eq 


0.538 


1997er 


0.466 


1997ct 


0.633 


1997CU 


0.592 


1997CX 


0.361 


1997ey 


0.575 


1997ez 


0.778 


1997fa 


0.498 



use the Kolmogorov-Smirnov (KS) statistic to test if a data 
set is consistent with the predictions of a specific model. 
Figure 1 (panels e-f) shows an example of this procedure. 
Bloom (2003) presents a similar treatment of the redshift 
distribution of gammar-ray bursts. 



3 APPLICATION TO AN OBSERVED SAMPLE 
OF HIGH-^ SUPERNOVAE lA 

We now choose a sample of SNe la whoso redshift distri- 
bution can be compared with predictions. Over the last 
few years, hundreds of high- 2; SNe la have been discovered, 
mainly by the SCP (Perlmutter et al. 1995) and the High-a 
Supernova Search Team (Schmidt ct al. 1998). These SNe 
were discovered during numerous observing runs, using dif- 
ferent telescopes and detectors, through various pass-bands 
and under non-uniform conditions. Given a complete de- 
scription of the observations, a combined analysis of all avail- 
able data sets should be feasible. At present, we limit our- 
selves to an analysis of a single data set obtained with a 
uniform observational setup, and for which most of the ob- 
servational details are available, in particular estimates of 
the survey efficiency r?. 

The SCP have published the results of their efficiency 
studies for several SN samples (Pain et al. 2002). Here, we 
study the deepest and largest of these data sets. This SN 
sample is denoted as set "D" in the analysis of Pain et al. 
(2002). Following is a brief outline of the search procedure 
used to discover these SNe. Deep (3600 s, total) 7?-band 
exposures were obtained with the BTC camera mounted 
on the 4m Blanco telescope at CTIO, on December 28-29, 
1997. These were compared with similar images obtained 
36-37 days earlier. Seventeen candidate SNe were discov- 
ered, with discovery magnitudes 21.6 < R < 24.5 and red- 
shifts 0.36 < z < 0.86. Twelve SNe were spectroscopically 
confirmed as SNe la. However, only eleven of these were re- 
tained in the "statistical sample" used by Pain et al. (2002) 
to calculate SN rates, and these constitute the sample we 
analyze, presented in Table 1. 

The SN sample was discovered in 11 different fields 
(Pain et al. 2002) observed under non-uniform atmospheric 
conditions, leading to survey efficiency curves r/(_R) which 
vaxy from field to field. We have averaged the efficiency 
curves measured for each field (kindly provided by R. Pain) , 
to get the mean efficiency of the survey, r]{R) (Figure 2). 




20 21 22 23 24 25 26 



R magnitude 

Figure 2. The mean efficieney of the SCP SN search described in 
§ 3 as a function of i?-band magnitude. Measurements of the de- 
tection efficiency in individual fields, kindly supplied by R. Pain, 
are given by -|- signs, while the open circles show the mean val- 
ues. The solid curve is a cubic spline fit to the mean points, and is 
used in our calculation. Note the large vertical dispersion resulting 
from non-uniform atmospheric conditions and from the variable 
fraction of each field that was not used due to the presence of 
bright stars or other defects. 



Since the efficiency enters our calculations in a linear fash- 
ion, calculating the SN redshift distribution using the mean 
curve is equivalent to performing the calculation for each 
field individually (with the appropriate efficiency curve) and 
summing the results. We neglect a possible redshift depen- 
dence of the efficiency due to the variation in the contrast be- 
tween SN light and underlying host galaxy emission within 
a constant aperture. This contrast is reduced with redshift 
oc DJ^^ (1-1-2)*, where Dl is the luminosity distance. For our 
chosen cosmology, this gives a factor of ~ 3 between z = 0.3 
and 2 = 1. Galaxy and SN K-corrections will offset some of 
this effect, and, from Fig. 2b of Pain et al. (2002) wc can 
estimate that the remaining correction to the SN detection 
efficiency will be of a few per cent at most. 

The search strategy employed by the observers in a cer- 
tain survey can also influence the predicted redshift distribu- 
tion of SNe la. In particular, the SCP survey whose data we 
use here requires that a candidate SN be brighter in a search 
image than in a reference image taken a few weeks earlier. 
This criterion was designed to discover only SNo near peak 
magnitude, and is driven by the SCP's scientific goals (Perl- 
mutter 1999, and references therein). The search criterion 
has implications for the calculation of the effective visibility 
times T{z). For example, low-2 SNe, that are much brighter 
than the survey's limiting magnitude Rum, could, in prin- 
ciple, contribute significantly to T (Eq. 8). However, during 
most of the time these objects spend above Rum, they are 
declining in magnitude, and therefore would have been ig- 
nored by the SCP. We therefore limit the calculation of T{z) 
to the period during which a SN is brighter in the second 
epoch image than it was 36 days before, when the reference 
images were obtained. This is achieved by multiplying the 
integrand in Eq. 8 by 
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"Madau" SFH 
(P=4. a=-2) 



redshift 




"Cowie-Steidel" SFH 
(P=1 , a=0) 



redshift 





Time deiay (x) [h ^ Gy] Time deiay (x) [h ^ Gy] Time deiay (x) [h ^ Gy] 



Figure 3. Probability of SN la time delay values, given the data, 
for particular SFH models. Assuming the SFH models shown in 
the upper panels, we can constrain the allowed values of r by the 
probability derived from the KS test (lower panels). Points below 
the upper and lower dotted lines are ruled out at 95% and 99% 
confidence, respectively. 



qtFH^^^^ - /l R(t)<R(t-36d) 

— \g R(t)>R(t-36d) 



(9) 



4 RESULTS 

4.1 Star formation history and the delay time of 
SNe la 

We have calculated the expected redshift distribution of 
SNe la for the above survey, for a range of characteris- 
tic delay times 0.3 < r < 3 Gyr, and SFH indices 
— 2 < Q < 2 and 1 < l3 < 4. Examination of the full range of 
models calculated shows that the data cannot constrain the 
model parameters individually, i.e., the characteristic delay 
time is degenerate with the SFH functional form. However, 
if we assume some prior, either on the time delay r or on 
the SFH, the data can be used to constrain the other pa- 
rameters. For instance, Figure 3 shows the range of allowed 
values for r assuming some of the SFH functions appearing 
in the recent literature. We can see that if the Madau et 
al. (1998) SFH is assumed, the delay time r is constrained 
to be T > 1.7(t > 0.7)h-^ Gyr at the 95%(99%) confi- 
dence level (CL). Rising SFHs, similar to those advocated 
by Lanzetta et al. (2002), are ruled out at the 95% CL, un- 
less T > 3h~^ Gyr. A gentler evolution of the SFH, e.g, as 
proposed by Cowie et al. (1999) and Steidel et al. (1999), 
places no constraints on t. Inspecting the actual distribu- 
tions predicted by the various models (e.g.. Fig. le and If), 
we find that invariably, the models that are ruled out by 
the data under-predict the number of SNe la at intermedi- 
ate redshifts (z ^ 0.5) relative to the number of SNe la at 
higher redshifts {z ~ 0.8). 

Figure 4 shows the constraints set on the SFH for given 
values of the delay time r. The first two panels show that, for 
short delay times (t < 0.5h~^ Gyr), the data constrain only 
the \ow-z index j3. For such short delays, the SFH at high z 
does not strongly affect (and is therefore not constrained 
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Figure 4. Constraints on SFH models for a given SN-Ia delay 
time. Each panel shows contours of equal probability, derived from 

the KS test for various combinations of the SFH indices a and f3, 
for the indicated SN la delay time, r, in units of Gyr. 



by) the SN redshift distribution we measure at 1. If we 
assume such short time delays, the SN data require a gentle 
decline in SFH from z = 1 to z = {f3 = 1), favoring the 

results of Cowie et al. (1999) and Sullivan et al. (2000a) over 
those of Lilly et al. (1996) or the recent findings of Hippelein 
et al. (2003). 

As we consider longer time delays (1 < t < 2h~^ Gyr), 
the data begin to constrain also the high-z evolution of the 
SFH, parameterised by the index a. Generally speaking, 
SFHs that decline above z = 1.2 (e.g., Madau et al. 1998) 
are favored over the strongly rising functions of Lanzetta et 
al. (2002). For instance, a model where the SFH is rising at 
all zas'^ ^ i'^ + zf (i-c, a = /3 = 2) is ruled out at the 95% 
CL for time delays r < 2h^^ Gyr. The constraints are weak- 
est when a long delay time, r = 3h~^ Gyr, is considered, in 
which case only SFH models combining a sharp increase at 
low redshift (/3 > 2) that continues to rise steeply at high 
redshifts (a = 2) are ruled out at the 95% CL, and are thus 
disfavored by the data regardless of r. In conclusion, even 
with a small sample of 11 SNe la, it appears that the SN la 
redshift distribution can be used to set interesting limits on 
the SFH and the characteristic delay time of SNe la. 

4.2 Caveats 

We now consider possible caveats that might relax the con- 
straints we have found. These may be due to some weakness 
in the model we have employed, or to effects related to the 
data set we have used. We discuss both possibilities below. 

Our model includes two main parameterised compo- 
nents the SFH and the SN la delay function. Our treatment 
of the SFH is fairly robust, since our parameterization covers 
the entire range of SFH functions reported in the literature. 
As for the SN la delay function, if SNe la are descended 
from a single or dominant type of progenitor, it will be rea- 
sonable to assume that there is a characteristic time delay 
between the formation of these systems and the SN explo- 
sion. Our parameterised models qualitatively reproduce the 
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Figure 5. As in Fig. 4, but using the DF time delay parameter- 
ization (§ 4.2). 

main features of most of the possible delay functions cal- 
culated for specific progenitor models by Ruiz-Lapuente & 
Canal (1998) and Yungelson & Livio (2000). A fastly rising 
and then monotonically declining delay function is also con- 
sistent with the observation that the rate of SNe la appears 
to be higher among young stellar populations than among 
older ones (e.g., Delia Valle & Livio 1994; Mannucci et al. 
2003), However, the functional dependence could be quite 
different. 

To further test the sensitivity of our results to the exact 
form of the delay function, we have redone our calculations 
using the alternative parameterization introduced by DF, 
who replace the exponential distribution of explosion times 
used by MDP by a discrete value of the delay r. The progen- 
itor mass then uniquely determines the delay time between 
the system formation and the SN la explosion. The results 
of the calculations axe shown in Figure 5. Comparing with 
Figure 4, we can see that all the main findings we obtained 
using the MDP delay functions are reproduced, and are in 
fact strengthened. As expected, because the delay time is no 
longer given by an extended function, the SN redshift dis- 
tribution between a ~ — 1 is influenced only by the low-z 
SFH index f5 for short delay times (r < \h~^ Gyr) and only 
by the high-z SFH index a for long delay times (t > 2h~^ 
Gyr). As before, for short delay times (r < 0.5/i.~^ Gyr) the 
data require a gently rising SFH {(3 = 1). For longer delay 
times (r > Ife"^ Gyr) we can rule out (at the 95% CL) SFH 
models with a sharp rise at high z (a = 2) while the models 
involving also a sharp rise also at low « (a = 2, /3 > 2) are 
ruled out at the 99% CL regardless of the time delay. 

Several studies (e.g., Hamuy et al. 1995; Hamuy et 
al. 1996; Howell 2001) have suggested that some "pecu- 
liar" SNe la (overluminous, SN 1991T-like, and underlumi- 
nous, SN 199 Ibg- like) may result from distinct subgroups 
of progenitors. It could therefore be argued that we should 
consider also models using two or more characteristic de- 
lay times. However, it has been shown that high- 2: SNe la, 
which dominate the sample we have considered, are at least 
as uniform, both spectroscopically and photometrically, as 
local "normal" SNe la (Riess et al. 1998; Perlmutter 1999, 



and references within). Indeed, Li et al. (2001) report that 
both overluminous and underluminous SNe have not been 
detected at high z so far. Overluminous SNe would be over- 
represented in high z samples. Their apparent absence prob- 
ably rules our the possibility that such events are common 
at high z. 

To check what effect an undetected population of un- 
derluminous events would have on our results, we have con- 
structed a set of SN light curves that have peak magni- 
tude and light curve shape similar to those of the prototyp- 
ical underluminous SN la, SN 1991bg. We then computed 
the effective visibility time using a mix of SNe la that in- 
cludes 30% 1991bg-like events, similar to the fraction found 
by Li et al. (2001) among local SNe la. Due to their rel- 
ative faintness, the redshift distribution of underluminous 
SNe la peaks around z ~ 0.4 (compared to z ^ 0.8 for nor- 
mal events. Fig Id). We find that this typically leads to an 
increase of ~ 15% to the total SN la distribution at this 
redshift. We noted above that the models calculated in sec- 
tion § 4.1, and which fail to fit the data, underpredict the 
number of SNe la at intermediate redshift {z ~ 0.4) relative 
to those at higher redshifts. Thus, the addition of SNe la at 
intermediate redshifts due to an underluminous SN popula- 
tion slightly weakens some the results reported above. For 
instance, assuming a Madau et al. (1998) SFH now requires 
T > 1.4(0.6)/i-^ Gyr at the 95%(99%) CL, while a Lanzetta 
et al. (2002) SFH can be accomodated if r > 2.5(1.2)/i-^ 
Gyr at the 95%(99%) CL. Other hmits are not relaxed - 
short (t < 0.5/i~^ Gyr) delay times require slowly varying 
SFH at 2: < 1.2 (/3 = 1), and steep, monotonically rising SFH 
(a = 2, /3 > 3) are ruled out at 95% confidence regardless of 
T. Thus, an undetected population of underluminous high-2 
SNe, if it exists, does not significantly alter our results. 

With regard to the data set we have analyzed, one may 
speculate that perhaps the observers excluded some SNe 
that should have been included in the sample, and may 
have thus affected the observed redshift distribution. How- 
ever, in their analysis. Pain et al. (2002) investigated all 
variable sources that could have been SNe la and which 
satisfied their search criteria. They determined whether or 
not each candidate was a SN la, based on spectroscopy and 
light curve shapes, leaving no unidentified sources. We have 
verified that the efficiency curves we have used are consis- 
tent with the selection criteria used by Pain et al. (2002) in 
compiling the SN sample shown in Table 1. 



4.3 Future Data Sets 

Finally, we consider potential applications of this analysis to 
future, larger, SN samples. We have created two simulated 
data sets. The first simulated data set consists of 100 SNe 
drawn from the expected distribution calculated for a survey 
with limiting magnitude R = 24, a "Madau" SFH (a = —2, 
(3 — 4), and a typical time delay of t = lh~^ Gyr. A second 
sample of 1000 SNe from a deeper survey to Rim = 25 
mag is produced using the same model parameters. We then 
treated each sample as real data and found constraints that 
could be drawn on the a — (3 — t parameter space. Samples 
similar to these synthetic sets are expected from ongoing or 
forthcoming SN search programs such as ESSENCE (Smith 
et al. 2002; Garnavich et al. 2002), SNLS (Pain et al. 2002), 
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Figure 6. As in Fig. 4, but for a simulated data sample of 100 
SNe, drawn from the distribution calculated for a = —2, 3 — 4 
and T = lh~^ Gyr and assuming a survey with limiting magni- 
tude R = 24. Note that the input model (marked with a filled 
circle) is recovered, and that significant areas of the parameter 
space are ruled out. However, individual limits either on the SFH 
indices {a and (3) or the delay time t cannot be drawn. 




Figure 7. Similar to Fig. 6, but using a simulated sample of 1000 
SNe to a limiting magnitude R = 25. Such a large data set could 
constrain both the SFH (/3 > 1 at the 95% CL) and the SN la 
delay time (t < 3h~^ Gyr) simultaneously. 



SNAP (Perlmutter et al. 2002) and the LSST (Strauss et al. 
2002). 

Figure 6 shows results of the first simulation. A data 
set of 100 SNe reaching E = 24 mag should be a powerful 
tool to test the SFH and SN la delay times. While the input 
model {a = —2, /? = 4, r = 1) is consistent with the data 
(P=0.49), significant portions of the parameter space are 
ruled out. However, the degeneracy between the SFH indices 
a and /3, and the SN la delay time t persists, and the data 
cannot constrain each parameter individually. 



Figure 7 shows the constraints that can be drawn 
from the larger, deeper sample. The "correct" model (a = 
—2,(5 = 4, T = 1) is recovered with a probability value of 
P = 0.79. Note that the degeneracy between the form of the 
SFH and the SN la delay time is now broken. Shallow SFHs 
at low z {(3 = 1) and long delay times (t = 3) are strongly 
ruled out by the data. Shorter delay times set strong con- 
straints on the allowed values of the SFH indices a and (3. 



5 DISCUSSION AND CONCLUSIONS 

We have calculated the expected redshift distribution of 
SNe la in a magnitude limited survey for a wide set of 
model parameters, and compared our calculations with ob- 
servations of high-z SNe by the SOP. We have shown that 
this approach is significantly more powerful than compar- 
ing expected and observed SN rates, because of the redshift 
binning involved in deriving SN la rates. We confirm and 
quantify the findings of Pain et al. (2002) and Tonry et al. 
(2003) which suggest that the redshift distribution of SNe la 
at 2 --^ — 1 disfavor a strong variation of the SN la rate 
during this period. Thus, if the SFH between z = and 1 is 
steeper than ^ oc (1 -|- z), the delay time of SNe la must be 
long (r > 0.5h~^ Gyr). We show that, even for longer delay 
times (r > Ih^^ Gyr), the data set interesting constraints 
on the SFH. Generally, strongly rising SFHs, similar to those 
advocated by Lanzetta et al. (2002), fit the data poorly, and 
SFH models which rise steeply at low redshifts (/? > 2) and 
continue to rise at high redshifts (a = 2) are ruled out by 
the data at the 95% CL, unless r > 3h^^ Gyr. Turning the 
argument around, given some determinations of the SFH 
(e.g., Madau et al. 1998), we can constrain the characteris- 
tic delay time of SNe la to be t > 1.7/i"^ Gyr (t > 0.7h~^ 
Gyr) at the 95% (99%) CL. 

As already noted, Yungelson & Livio (2000) have stud- 
ied the delay functions of different SNe la progenitor sys- 
tems. In their double-degenerate (DD) model, two WDs 
merge, with the resulting WD reaching (or surpassing) the 
Chandrasekhar mass. Single degenerate models include ac- 
cretion of He-rich material from a non-degenerate compan- 
ion, leading to helium ignition on the WD surface and an 
edge- lit detonation (He-ELD). Alternatively, H-rich material 
from a main sequence companion is accreted and processed 
on the WD surface, and results either in an edge-lit deto- 
nation due to accumulated helium (MS-ELD) or to central 
carbon ignition as the WD reaches the Chandrasekhar mass 
(MS-CH). 

The delay functions calculated for these models are 
more complex than the parameterised MDP forms we have 
used, but generic similarities exist. In particular, both DD 
and He-ELD curves can be characterised by a fast onset of 
SNe la, some 3 x 10^ years after star formation, followed by 
gradual increase and a decline that terminates in a sharp, 
exponential cutoff after ~ 1 Gyr for the He-ELD and ~ 11 
Gyr for the DD models. This is qualitatively similar to the 
MDP exponential formulation we use. The MS-CH model 
can be approximated by a DF-like model with a delay time 
of ~ 1 Gyr, so can be constrained by our calculation using 
this parameterization (§ 4.2, Figure 5). The MS-ELD is more 
complex, with a delayed onset of SNe la, about 0.3 Gyr af- 
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ter star formation, followed by an approximately power-law 
decline lasting some 11 Gyr. 

If we consider those SFH determinations requiring long 
time delays in our models (e.g., Madau et al. 1998; Lanzetta 
et al. 2002), then the data disfavor the MS-CH and He-ELD 
progenitor models. We have shown that larger samples of 
SNe la, that will be obtained by future SN search programs, 
some of which are already in progress, could be used to ob- 
tain even stronger constraints on the SFH and SN la models. 

In a companion paper (Maoz & Gal- Yam 2003), we 
study the implications of the measured SN la rates in 2; ~ 1 
galaxy clusters on the nature of SN la progenitors and the 
source of iron in clusters. We demonstrate there that, if 
SNe la have produced most of the iron in galaxy clusters, 
and the stars in clusters formed at 2; ~ 2, then the SN la 
delay time must be lower than 2 Gyr. Thus, if both condi- 
tions were met, then the Lanzetta et al. (2002) SFH would 
be ruled out by the data presented here. 

Finally, it is interesting to note here that most previ- 
ous authors have considered core-collapse SNe, whose rates 
closely track the SFH, as promising indicators for SFH mear- 
surement (but see DF 1999 for an alternative view). In prin- 
ciple, our methods could easily bo applied also to liigli-z 
core-collapse SN samples, if and when they become avail- 
able. However, recent work by Mannucci et al. (2003) has 
shown that the majority of core-collapse events in star- 
forming galaxies suffers from strong dust extinction. SFH 
measurements using core-collapse SNe thus share the same 
problems that make UV-determined SFHs so widely debated 
- the strong dependence of the derived SFH on the little 
known amount and properties of dust at high redshifts. In 
contrast, most SNe la probably occur in relatively dust-free 
environments. Thus, measurements of the SFH using SNe la, 
through methods similar to the ones we have discussed here, 
may provide an attractive alternative to SFHs based on UV 
and emission line fluxes. 
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